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We investigate the specific heat and magnetisation of a 
ferrimagnet with gS = 1 and S — | spins in a finite magnetic 
field using the transfer matrix DMRG down to T = 0.025J. 
Ferromagnetic gapless and antiferromagnetic gapped excita- 
tions for H = lead to rich thermodynamics for H > 0. 
While the specific heat is characterized by a generic dou- 
ble peak structure, magnetisation reveals two critical fields, 
H c i — 1.76(1) and H C 2 = 3.00(1) with square-root behaviour 
in the T = magnetisation. Simple analytical arguments 
allow to understand these experimentally accessible findings. 



In recent years, one-dimensional ferrimagnets with al- 
ternating|-|S»ins gS and S have attracted considerable 
attentionEnll The Lieb-Mattis theorem reveals the char- 
acteristic feature of ferrimagnets: the ground state of a 
ferrimagnet with N elementary cells of two spins is a 
macroscopic spin of length NS(g — 1). While the inter- 
action is antiferromagnetic, the ground state resembles 
that of a ferromagnet with Neel-like alignment of big and 
small spins. Due to quantum fluctuations, the classical 
Neel state is not exact, but the macroscopic magnetisa- 
tion of the ground state makes spin wave theory applica- 
ble, even in one dimension. 

Spin wave theoryo yields two types of excitations: 
starting from magnetisation NS(g — 1), there are fer- 
romagnetic (FM) gapless excitations to states with mag- 
netisation NS(g — 1) — 1 (Goldstone modes), and anti- 
ferromagnetic (AFM) gapped excitations to states with 
magnetisation NS(g — 1) + 1. 

These two excitation types lead to a crossover in the 
beJiaviour of the specific heat C v and the susceptibility 
xO'El. C v shows an AFM mean field peak at intermediate 
temperatures and a FM C v cx VT behaviour for T — > 0. 
X oc T _1 for T — > oo, but shows a FM T~ 2 divergence at 
T ^ 0. 

At finite field all dispersion relations are shifted and 
the ground state degeneracy is lifted, such that experi- 
mentally easily observable changes in the thermodynamic 
quantities should occur. In this paper we discuss the ex- 
pected observations and provide very precise quantitative 
results obtained mainly by the transfer matrix Density 
Matrix Renormalization Group (DMRG)Lm, hoping to 
stimulate further experimental investigation. 

We consider the generic (1, h) ferrimagnet with gS = 1 
and S = h. For this case, the AFM gap in zero field is 
numerically found! to be A AFM (0) = 1.759 J. 

Analytical results. The introduction of a field term 



H J2i Si leaves all eigenstates invariant, while shifting the 
eigenenergies by HS% ot , where 5f ot is the total magneti- 
sation of the eigenstate. 

Linear spin wave analysis based on the classical fer- 
rimagnetic ground state gives the following two disper- 
sion "branches" (all energies, temperatures and fields are 
measured in units of J = 1): 



(q)= v /5-4cosg±l )/2tH 



(1) 



Spin wave theory gives eigenstates of S% ot ; thus the 
only effect of an external field on a one-particle exci- 
tation carrying magnetisation ±1 is the Zeeman term 
introducing a gap A FM (H) = H to the FM excita- 
tions cu~ , while the AFM branch lo + is reduced in gap: 
A AFM {H) = A AFM (0)-H. 

In linear spin wave theory, m-particle excitations car- 
rying magnetisation ±m are simply linear superpositions 
of m one-particle excitations. As the field acts only 
through the Zeeman term, it is sufficient in this approxi- 
mation to consider a one-particle excitation to study field 
effects. 

At H = A AFM (0), AFM excitations become gapless, 
allowing spin-up flips at no energy cost, such that the 
simple result M(T = Q) = 0.5 breaks down. Using the 
numerical value for the gap, A = 1.759, rather than the 
spin wave result, we conclude that there is a first critical 
field H c i = A(H = 0) = 1.759, above which M(T = 0) > 
0.5. 

Second, for a field H = A(H = 0)/2 (i.e. H cl /2, 
low temperature magnetisation M(T) should be con- 
stant with T, while it is monotonically decreasing with 
T for H < A/2, and monotonically increasing with T for 
H > A/2. This is because the dispersion relations of the 
harmonic spin waves increasing and decreasing magneti- 
sation become identical at this field value. 

Linear spin wave analysis based on the fully polarized 
state yields a dispersion relation 



to ± {q)=H 



^±^5 + Acosq = H-LO F± (q), (2) 



where uj 1 (q) is the two-branched exact FM one-particle 
excitation: both excitation branches reduce magnetisa- 
tion. Below a second critical field H c2 = 3, one dispersion 
branch acquires negative energy, such that a breakdown 
of the simple result M(T = 0) = 3/2 is predicted for 
H < H c2 . 

Let us now consider the specific heat. C V (T) should 
acquire a generic double-peak structure for low fields < 
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H < H c i, because two gapped excitation modes exist: 
Both gapped antiferromagnets as well as ferromagnets in 
external fields (see Fig. 0) show exponential activation 
of C v with a pronounced peak, whose position is related 
to the gap or field size. Peak positions should thus shift 
with H to higher (FM contribution) and lower (AFM 
contribution) temperatures because of the Zeeman term. 
At the critical fields, where the gapless excitation has the 
form w oc g 2 , C v cx \/T is expected for T — > 0. 

Another analysis for low fields can be obtained from 
studying the decoupled-dimcr limit, useful in the H = 
caseEl: every second interaction is switched off, yielding 
a flat dispersion, but the FM and AFM elementary ex- 
citations are still gapless and gapped, respectively. This 
analysis gives a double-peak structure for C v . With in- 
creasing field, the double-peak structure is smeared out 
and the two peaks merge into one. With further increase 
of the field, the system changes discontinuously into the 
fully polarized state. 

The existence of two critical fields H c i (up to which 
the ground state magnetisation persists) and H C 2 (which 
marks the beginning full polarisation) in the ferrimag- 
net is analogous to the finite field behaviour of the 5 = 1 
(Haldane) chain. In both systems the characteristic prop- 
erty of gapped excitations existing for both < H < H c \ 
and for H > H C 2 breaks down at the critical fields; in 
spin wave theory this is indicated by an instability of 
an arbitrary number of spin waves (spin wave condensa- 
tion). As for the S = 1 chain (Ref. 0) we expect the oc- 
curence of a critical (Luttinger liquid) phase in the inter- 
mediate regime H c \ < H < H C 2- This phase should re- 
semble a (critical) anisotropic Heisenberg chain with the 
anisotropy determined by the magnetic field; power law 
correlation functions as well as a linear behaviour of the 
specific heat C v — jT (with 7 related to H) are expected. 
The ferrimagnetic chain, however, having two spins per 
unit cell, is richer than the 5=1 chain and the inter- 
play of two elementary excitations leads to new crossover 
phenomena between high and low field behaviour in both 
specific heat and magnetisation. 

The spin wave approach thus results in an interesting 
qualitative picture, but evidently it suffers from a num- 
ber of deficiencies: The results for the non-Zeeman part 
of the excitation energies from the ferrimagnetic ground 
state is only approximate, the stiffness of the FM exci- 
tations is overestimated and the gap energy of the AFM 
excitations is underestimated (Ref. ||), for quantitative 
estimates the numerical values should be used. Spin wave 
theory does not yield a peak for the specific heat as the 
harmonic approximation implies dC v /dT > for an arbi- 
trary density of states; however, for small fields, it gives 
an unusually high C v at very low T, consistent with a 
second low-T peak. 

The predictions of spin wave theory for the critical 
fields as described above, further rest on the assumption 
that multi particle excitations do not become unstable 
before the one particle excitation. If this is true, as is 
plausible and generally asumed, the value H c2 = 3 for the 



upper critical field is exact (since the corresponding one 
particle excitation energy is exact); this appears likely 
although a formal proof is not available. 

Transfer Matrix DMRG Results. To obtain quantita- 
tively reliable values for the specific heat and the mag^ 
netisation, the recently proposed transfer matrix DMRGlI 
was applied to the problem. Essentially a decimation pro- 
cedure applied to a quantum transfer matrix, it maintains 
the advantages of the quantum transfer matrix method 
such as working in the thermodynamical limit, and al- 
lows the evaluation of very large Trotter numbers, giving 
reliable access to thermodynamics at very low tempera- 
tures. In the problem under study, most interesting ob- 
servations can be made at very low temperature; in the 
critical region H c \ < H < H C 2 the absence of finite size 
effects is of advantage. 

The (controlled) approximation of the DMRG rests 
on keeping a reduced state space. We find that keep- 
ing M — 80 states yields almost exact results (note that 
we overestimate the specific heat by several percent at 
higher temperatures, due to coarse DMRG sampling at 
high T: sampling temperatures are spaced by A/3 = 0.2), 
excluding the case H — 0, where the very high FM de- 
generacy makes it necessary to push the method much 
furtheiO. 

In the specific heat per elementary cell, a two- 
peak structure evolves for small fields, with the high- 
temperature AFM peak moving to smaller temperatures 
and the low-temperature FM peak moving to higher 
temperatures (Figure |l|). Using for a first rough es- 
timate C v cx x 2 cosh - x, x — JA/2T, leading to 
A w 2AT pea k, one finds for the high-T peak a zero-field 
gap A cv 1.70J in reasonable agreement with the nu- 
merically observed gap (A sw 1.759J); it shifts linearly 
and almost proportionally with H to lower T. As Fig. 
H shows, the low-T peak is almost identical with that of 
a 5 = i ferromagnet; the "remaining" specific heat can 
be attributed to AFM excitations (this naive comparison 
is of course not perfect, but renormalization effects are 
surprisingly weak). For H > 0.2, the two peaks merge 
into one; a shoulder for H = 1.6 indicates the emergence 
of a second peak (Figure ||). For H > H cl a two peak 
structure emerges again (Figure ||). In the intermediate 
regime H c \ < H < H C 2 the expected linear behaviour 
of the low temperature specific heat is clearly observed 
for fields H — 2.0, 2.4. For magnetic fields close to the 
critical values our results seem to indicate a behaviour 
C v cx VT for T — > 0. 

Considering the magnetisation per elementary cell, one 
finds for H < H cX that M = 1/2 at T = 0. How this 
magnetisation is reached for a fixed field H, is due to 
a competing effect of magnetising interactions and ther- 
mal fluctuations. For high temperature, thermal fluctu- 
ations will always suppress magnetisation. For fields be- 
low H = 0.88, demagnetising FM fluctuations cost less 
in energy and also suppress magnetisation. For higher 
fields, magnetising AFM fluctuations will increase mag- 
netisation. One observes therefore an intermediate mag- 
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netisation peak before thermal fluctuations again sup- 
press magnetisation (Figure ||). We observe an almost 
flat magnetisation curve up to T « 0.2 at H 0.88. in 



exce llent accordance with the prediction (Figure 



We find the lower critical field (Figure |7]) at H c \ = 
1.76(1), the ferrimagnetic magnetisation at T — is bro- 
ken up at H = Hoi, with M (H) - M (H cl ) cx y/H - H cl 
(cf. Refs. ^ for equivalent observations in Haldane mag- 
nets). For T -> 0, we find that (M(H) - M(H cl )) 2 
becomes linear in H close to H cl . As the system is 
critical at this field, working with the transfer matrix 
DMRG, which is in the thermodynamic limit, is prefer- 
able to conventional DMRG, which would suffer frop, 
finite-size effects and possible metastable trappingst2l 
due to extensive level crossings at H c \ and immediately 
above. Numerically, we locate the upper critical field at 
H C 2 = 3.00(1), in excellent agreement with the analytical 
result. Again, a square root behaviour in the approach 
to the maximum magnetisation can be observed. 

Conclusion. We have shown that the properties of a 
ferrimagnet in a weak external field can be understood 
qualitatively and in some cases even quantitatively in 
an extremely simple picture. The specific heat clearly 
reflects the dual structure of the excitations of the ferri- 
magnet for all fields. The magnetisation is dominated by 
the two critical fields, which can be understood as pure 
Zeeman effects. The Luttinger-like critical phase in be- 
tween, with square root divergences in the magnetisation, 
needs consideration of many-body effects for understand- 
ing. The presented finite-temperature properties, in par- 
ticular the low-field behavior of C v should be accessible 
to experimental verification. 
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FIG. 1. C v vs. T for fields H < 0.4. 
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FIG. 2. C v at H = 0.05 for a ferrimagnet (dotted, DMRG) 
and a S = i ferromagnet (solid, DMRG). To a very good 
approximation, the difference in specific heat (dashed) can be 
identified as AFM contribution. 
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FIG. 3. C v vs. T for H close to H cl = 1.76(1). 
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FIG. 4. a vs. T for fields H > H cl = 1.76(1). 
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FIG. 5. Magnetisation vs. T for fields H = 0.1, 0.4, 0.8, 
1.2, 1.6, 1.75, 1.9, 2.4, 2.8, 2.9, 3.0, 3.1, 3.2 (from the lowest 
to the highest curve) . 
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FIG. 6. T-independence of the low-T magnetisation for 
magnetic fields close to half the critical field H c i — 1.76(1). 
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FIG. 7. Magnetisation vs. H for T = 0.025, 0.05 (DMRG), 
T — 0.08, 0.1, 0.15 (quantum MC, to show consistency of 
methods). Curves approach singular behaviour for T — > 0. 
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